New limit on the K\ —> n®vv decay rate 
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The first dedicated experiment for the rare kaon decay K\ — > i^vv has been performed by the 
E391a collaboration at the KEK 12-GeV proton synchrotron. A new upper limit of 2.1 x 10 -7 at 
the 90 % confidence level was set for the branching ratio of the decay K\ — > i^Pvv using about 10 % 
of the data collected during the first period of data taking. 

PACS numbers: 13.20.Eb, 11.30.Er, 12.15.Hh 



The rare decay if£ — > ifivv is a Flavor Changing 
Neutral Current (FCNC) process from strange to down 
quarks and is caused by direct CP violation [J, Q. 
The theoretical uncertainty in the branching ratio is 
only 1-2%, while the branching ratio is predicted to be 
(2.8 ± 0.4) x 10 -11 3] based on the currently known pa- 
rameters from other experiments. The decay is consid- 
ered an ideal process to study the quark flavor physics 
and a critical test of the standard model as well as a 
search for new physics beyond it j4j . The present experi- 
mental limit is 5.9xl0 -7 at the 90 % confidence level Q; 
the Dalitz decay mode 7r° — > e + e~7 for the final state of 
K\ — » ifivv was used in the search. 

The E391a experiment at the KEK 12-GeV proton syn- 
chrotron was proposed to be the first dedicated experi- 
ment for the K\ — > ifivv decay and aimed to improve the 
experimental sensitivity by orders of magnitude, and to 
verify the experimental method for the next higher sen- 
sitivity experiment 6]. The E391a experiment had three 
data-taking runs in 2004 and 2005. The first data taking 
(Run-1) was performed from February to June 2004 [g. 
In this Letter, we report results obtained from an analysis 
of about 10 % of the data collected in Run-1. 

The signature of the K\ — > tt°vv decay is 27 + noth- 



ing in the final state. The energies and hit positions 
of two photons were measured with an electromagnetic 
calorimeter. The 7r° vertex, Z vtx , and its transverse mo- 
mentum with respect to the beam axis, Px, were mea- 
sured assuming that two photons were produced in a 
7T° — > 77 decay on the beam axis. One of the crucial tools 
for this detection method is a neutral beam with small 
diameter (called "pencil beam" ) in order to minimize un- 
certainty in the Z vtx and Pt measurements. Good beam 
collimation with little halo is also important to minimize 
7T° production via interactions with detector material. 
The second tool is a hermetic veto system covering the 
decay region to reject background decay modes with ad- 
ditional charged particles and photons. High efficiency 
can be obtained by lowering the energy thresholds |£j . 
The third tool is a decay region with a high vacuum to 
reduce 7r° produced by neutrons in the beam interacting 
with the residual gas. 

Neutral kaons were produced by striking protons from 
the 12 GeV proton synchrotron on a 60 mm long plat- 
inum target. The neutral beam was collimated to a circu- 
lar shape having 2 mrad of a half cone angle at an angle 
of 4° with respect to the primary proton line [To| . 

Figure □ shows the E391a detector. The "0 m" in the 
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figure corresponds to 11 m downstream from the target. 
Detectors were cylindrically assembled along the beam 
axis, and most of them were installed in a large vacuum 
vessel to eliminate dead material in front of the detectors. 
The fiducial region for decays was from 300 cm to 500 
cm. 

The electromagnetic calorimeter was placed at the 
downstream end of the decay region, and was 1.9 m in 
diameter with a 12 x 12 cm 2 beam hole at the center. 
It consisted of 576 blocks of undoped Csl crystal, with 
a size of 7 x 7 x 30 cm 3 (=16Xo) except for the in- 
ner 24 crystals that were 5 x 5 x 50 cm 3 (=27X ) [ll|. 
The energy resolution was cje/ E ~ 1%/^/E 01%, where 
E is in GeV, as measured with 25 Csl blocks and a 
positron beam. The average position resolution was 5 
mm for photons. A group of plastic scintillation coun- 
ters (CV) were placed in front of the Csl calorimeter to 
veto decays involving charged particles. The decay re- 
gion was surrounded by two large lead-scintillator sand- 
wich counters, MB and FB. Total thickness of MB (FB) 
was 13.5 Xo (17.2 Xo). Their signal was read through 
wave-length-shifting fibers by phototubes (PMTs) with 
high quantum efficiency . The light yield of MB and 
FB was monitored with an LED calibration system. The 
gain shift between the on-beam period and the off-beam 
period was less than 1 % for both MB and FB. Multiple 
collar-shaped veto counters, CC02 - CC07, were placed 
along the beam axis to detect photons going through the 
beam hole. CC02 was a shashlik type lead-scintillator 
sandwich counter and CC03 was a tungsten-scintillator 
counter with the layers set parallel to the beam. CC04 
and CC05 were lead-scintillator counters with the layers 
set perpendicular to the beam. CC06 and CC07 con- 
sisted of SF-5N lead-glass blocks. Back Anti (BA) was 
located at the end of the beam in order to detect pho- 
tons going through the beam hole without being detected 
by other detectors. BA consisted of lead-scintillator lay- 
ers and quartz layers, and had 14 Xo in total. In front 
of BA, a set of 1 mm thick plastic-scintillators, BHCV, 
was placed to detect charged particles going through the 
beam hole. 



The vacuum tank was divided internally into two re- 
gions by a 20 mg/cm 2 thick sheet called "membrane" 
to protect the high vacuum region from out-gassing from 
detector components. The pressure of the region with de- 
tector components was less than 1 Pa, and the pressure 
in the decay region was 1 x 10 -5 Pa. 

For the Csl calorimeter and all the photon veto de- 
tectors, we continually calibrated the energy scale factor 
using cosmic ray muons and minimum ionization parti- 
cles in the beam during the operation. We also studied 
the energy scale factor of the Csl calorimeter using the 
special data, in which an aluminum(Al) target with a 
thickness of 5 mm was inserted in the beam at z = 280.5 
cm (downstream of CC02) in order to produce 7r° by 
neutrons in the beam. We adjusted the gain of each Csl 
crystal so that the invariant mass of two photons agree 
with the 7T° mass. This data was also used to check our 
event reconstruction performance. The more detailed de- 
scription of the calibration can be found in 

The trigger was designed to accept K\ — > 7r ^, 
K° L 7r°7r° and K° L 7r 7r°7r° decays. The latter two 
decays were used to study the detector response, and the 
K\ — > 7r°7r° decays were used to measure the number 
of K L decays. For triggering purpose, we grouped eight 
neighboring Csl crystals and defined 72 regions in total. 
The trigger required that there were two or more such 
regions with > 60 MeV energy deposit in each. The trig- 
ger also required no energy deposit in CV and several 
photon veto detectors. For example, the veto threshold 
for the total energy deposit in MB was 15 MeV. The trig- 
ger rate was 800 events per 2 second beam delivery with 
a typical intensity of 2.5 x 10 12 protons on the target. 
The live time ratio was 78 % with a network distributed 
data acquisition system with multiple CPUs. The elec- 
tronics and data acquisition system is briefly described 
elsewhere [?1 Hlj. 

In the offline analysis, we first looked for photons in 
the Csl calorimeter. Each cluster of energy deposits was 
required to have the transverse shower shape consistent 
with a single electromagnetic shower. The effective en- 
ergy threshold of each cluster was 10 MeV. We assumed 
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that the two photons came from a 7r° decay, and recon- 
structed Z vtx requiring the two photon invariant mass to 
have the 7r° mass. 

We selected events with exactly two photons hitting 
the Csl calorimeter and applied selection criteria (cuts) 
to suppress background events. In Run-1, the down- 
stream membrane was partially hanging in the beam by 
error, at z = 550 cm. This produced a large number 
of background events because neutrons in the beam core 
struck the membrane and produced secondary 7r 0, s. If 
multiple 7r 0, s were produced at the membrane, and two 
photons from different 7r 0, s were detected ("core neutron 
multi-7r°" event), it became a serious background event 
because we were not able to reconstruct Z vtx correctly 
and these events were distributed in the fiducial region. 
On the other hand, these events had extra photons in 
the final state, and thus can be suppressed by detecting 
those extra photons. 

In order to suppress events involving extra photons, 
we required energy deposit in each photon veto detec- 
tor to be less than the threshold listed in Table [I] The 
rejection power of the photon veto was evaluated with 
four-photon event samples from K\ — > 7r°7r° — > 47 and 
if£ — > 7r 7r°7r° — > 67 with two missing photons. Fig- 
ure El shows the invariant mass of four photons, M~4 7 , 
after applying all the cuts on the photon veto detectors. 
With all the photon veto cuts, the ratio of the number of 
Kl 7r°7r° events in 0.45 < M 47 (GeV/c 2 ) < 0.55 to the 
number of K\ -> 7r 7r°7r° events in M 47 (GeV/c 2 ) < 0.45 
improved by a factor of 1 1 . This improvement was consis- 
tent with the expectation of GEANT-3 based Monte 
Carlo simulation (MC) within 18 %. 



TABLE I: List of the thresholds applied to the photon veto 
detectors. Eq is the total light yield in the BA quartz layers, 
and Es is the total energy deposit in the BA scintillator layers. 
The signal efficiency for a cut A, sa, is the ratio of the number 
of events with all cuts to the number of events with all cuts 
except for the cut A. We estimated sa using MC K\ — > i^vv 
events except for the cut on BA. For BA, we first evaluated ea 
for K\ — > 7r°7r° and K\ — > 7r 7r°7r° decays using real data, which 
were 0.638=b0.022( stat- ) and 0.658=b0.022( stat- ) , respectively, and 
then assigned the average as sa for K\ — > i^vv signal. 



Detector threshold sa Detector threshold 


SA 


CC02 


4 MeV 1.0 CC06 5 MeV 


0.98 


CC03 


1.5 MeV 0.98 CC07 50 MeV 


0.99 


CC04 


3 MeV 0.98 FB 2 MeV 


0.91 




Detector 


threshold 


SA 


Csl 


3 MeV for the Csl crystals which do 
not belong to the photon clusters 


0.78 


MB 


1 MeV for the inner modules, and 
0.5 MeV for the outer modules 


0.60 


BA 


0.5 MIP for E Q , and 

Eq/Es > 10 


0.65 




FIG. 2: Distribution of the invariant mass of four photons, 
M4 7 , with the cuts on the photon veto detectors and the 
shower shape of photons in the Csl calorimeter. The dots 
show data, the open solid histogram shows total MC, the 
closed solid histogram shows K L — > MC and the hatched 
histogram shows K\ 



From the MC study, we found that the "core neu- 
tron mufti- 7T°" events had low Pt and were populated 
at the downstream Z vtx region. In order to minimize the 
number of such background events, we used this charac- 
teristic and required a parameter, a = Pt (GeV/c) — 
8.0 x 10- 4 (GeV/c-cm) x Z vtx (cm) (Z vtx < 525cm), to 
be larger than -0.225 GeV/c. Another cut, P T > 0.12 
GeV/c, was applied to suppress K\ — > 77 background 
and A — > ir° n background, whose maximum Pt is 0.109 
GeV/c. The upper boundary on Pt was determined to 
be Pt < 0.24 GeV/c from the kinematical limit of the 
K\ — > ifivv decay (P m ax = 0.231 GeV/c), allowing for 
the smearing effect due to detector resolutions. 

After applying all the selection cuts, we estimated the 
number of remaining background events in the eight Pt~ 
Zvtx regions with the signal regions (c) and (d) as shown 
in Fig. [3 Except for the regions (a), (c) and (g), the dom- 
inant background source was the "core neutron multi-7r°" 
event. We evaluated the number of "core neutron multi- 
7T°" events using a relational expression with two inde- 



pendent selection cuts: N^ g = N f x (cut-1 rejection) x 
(cut- 2 rejection), where "cut-1" is a set of cuts on CV, 
MB, CC03, CC04, CC06 and CC07, "cut-2" is a set of 
cuts on the cluster energy and the cluster hit position, 
and N' is the number of events with all the selection cuts 
except for the cut-1 and the cut-2. We checked that (i) 
(97 ±3) % of the N' was the "core neutron multi-7r°" 
events even without the cut-1 and cut-2, and (ii) the cut- 
1 and the cut-2 were independent of each other 15]. The 
number of "core neutron multi-7r°" background events 
was O.Oio.o m re gi° n ( c ) an d 1.5 ± 0.7 in region (d). The 
number of background events caused by the halo neu- 
trons interacting with the detector material (CC02) and 
producing one or more 7r° was 0.9 ± 0.2 in region (a) 
and 0.04 ± 0.04 in region (c). The background events 
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FIG. 3: Zvtx versus Pt with all the event selection cuts. The 
number of observed (total expected background) events are 
shown. The expected number of background events was con- 
sistent with the observed number of events for all the regions. 



caused by the core neutrons interacting with the mem- 
brane and producing z/'s ("core neutron rf events) was 
reconstructed around region (c) because the 7r° mass was 
assumed. For all the events, we re-calculated the decay 
vertex assuming r] mass (Z v ) and then rejected events 
around the membrane in the beam, 525 < Z v (cm) < 575. 
The remaining number of "core neutron 77" events was 
0.4 ± 0.2, which was the largest component in region 
(c). The background events from K\ — > 7r°7r° with two 
missing photons was evaluated with MC. The number of 
K\ — > 7r°7r° background events in the signal region was 
0.04 ± 0.03, where the error includes the MC statistics 
and the systematic uncertainties, of which the dominant 
source was the mismatch between data and MC in the 
transverse shower shape of photon in the Csl calorime- 
ter. Moreover, the K\ — > 7r°7r° background events were 
the largest component in region (g). The total number 
of background events in the signal region was estimated 
to be 0.4^2 in region (c) and 1.5 ± 0.7 in region (d). 

We estimated the acceptance of K\ — > ifivv decay 
to be (0.657 ± 0.016) x 10 -2 based on cut efficiencies 
evaluated with the real data and MC study. The main 
components of the acceptance loss were the cuts on MB 
and BA photon veto detectors. In order to estimate the 
number of K\ decays in this search, we analyzed K\ — » 



7r u 7r u decays. The invariant mass and the reconstructed 
decay vertex for — > 7r°7r° are shown in Fig. In the 
signal region: 0.47 < M 47 (GeV/c 2 ) < 0.53, 




decay vertex (cm) 



FIG. 4: Distribution of the invariant mass (left) and the decay 
vertex (right) for the K\ — > 7r°7r° decays. In the top plot, the 
dots show the data and the histogram shows the MC. The 
bottom plot shows the ratio of the data to the MC. 



and 300 < Z vtx (cm) < 500, there were 2081 K% -> 7r°7r° 
events after subtracting 30 K\ — > 7r 7r°7r° background 
events. Based on the MC study, we estimated that the 
acceptance of K\ — > 7r°7r° decay was 1.41 x 10 -3 . 

The different final states between the signal and nor- 
malization modes caused systematic uncertainties in the 
single event sensitivity. We assigned the total systematic 
uncertainty in the single event sensitivity to be ± 7.0 %. 
The large sources of systematic uncertainty came from 
the mismatch between data and MC in the transverse 
shower shape of the photon (4 %) and the energy distri- 
bution in MB (4.2 %). 



With the K° L 



A 



7r°7r° branching ratio, (8.83 ± 0.08) x 
10 -4 14], we estimated the number of K\ decays to be 
(1.67±0.04( stat> )) x 10 9 . The single event sensitivity was 
(9.11±0.20 (stat .) ±0.64 (syst .)) x 10~ 8 . Since we observed 
no events in the signal region, we set a new upper limit on 
the branching ratio of K\ — » ifivv to be < 2.1 x 10 _T at 
the 90 % confidence level based on the Poisson statistics. 
This represents an improvement of a factor of 2.8 over 
the current limit p. 
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